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ABSTRACT We deduce the dynamic frequency-domain-lattice Green's function of a linear chain with 
properties (masses and next-neighbor spring constants) of exponential spatial dependence. We analyze the 
system as discrete chain as well as the continuous limiting case which represents an elastic ID exponentially 
graded material. The discrete model yields closed form expressions for the N x N Green's function for an 
arbitrary number N = 2,..,oo of particles of the chain. Utilizing this Green's function yields an explicit 
\^ | expression for the vibrational mode density. Despite of its simplicity the model reflects some characteristics 
if} ■ of the dynamics of a ID exponentially graded elastic material. As a special case the well-known expressions 
for the Green's function and oscillator density of the homogeneous linear chain are contained in the model. 
The width of the frequency band is determined by the grading parameter which characterizes the exponential 
| spatial dependence of the properties. In the limiting case of large grading parameter, the frequency band 
is localized around a single finite frequency where the band width tends to zero inversely with the grading 
parameter. In the continuum limit the discrete Green's function recovers the Green's function of the continuous 
equation of motion which takes in the time domain the form of a Klein-Gordon equation. 
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1 INTRODUCTION 

There is an increasing interest in many dynamic applications to control the frequency band of vibrational 
modes. Extensive models dealing with how point defects in crystalline lattices affect the density of vibrational 
modes by using lattice Green's functions [1]. As a rule in these models the defects are conceived as sufficiently 
"small" perturbations of a regular lattice and treated in the framework of a perturbation calculus where the 
lattice is considered in the harmonic approximation [2, 3, 4]. All these models have in common that they 
consider the effect of localized defects on the vibrational spectrum resulting in the appearance of "soft modes" 
which can be verified by their effects on the low temperature behavior of the specific heat [4, 5]. 
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However, little is known about the vibrational mode density when the lattice has "strongly" spatially 
varying properties which cannot be conceived as "small perturbations" with respect to an homogeneous 
reference lattice. In recent decades the study of spatially inhomogeneous materials has become of increasing 
interest, in particular due to their increasing importance in new applications in so called Smart Materials and 
as functionally graded materials (FGM) . However their technological exploitation still is limited in view of the 
principal difficulty to set up analytic continuum models. A principal reason for this is the fact that even for 
their simplest variants, the FGM with exponential properties, continuous Green's functions are not available 
in explicit forms. This is true for static Green's functions [6] and more than ever for the dynamic framework 
such as as for problems of heat conduction [7] in two and three dimensions. 

The present paper aims at establishing a simple discrete one-dimensional linear chain model for exponen- 
tially graded material. The exponential spatial dependence of the properties, i.e. of the masses and spring 
constants is assumed such that their ratios remain spatially independent. The frequency-domain dynamic 
lattice Green's function of this ID system is determined in closed form for the entire range of frequencies, 
that is within the band of eigenfrequencies and outside the band of frequencies. The expression deduced 
for the N x N Green's function matrix is explicit for an arbitrary number N > 1 of particles in the chain. 
Utilizing this Green's function gives also an exact expression the density of oscillation modes. It turns out 
that variation of the frequency band as well as the lowest and highest frequency can be controlled by the 
grading parameter. In spite of the simplicity of the model, the latter provides some of the essential features 
of one dimensional FGM. 



2 The Discrete Model 

We analyze a linear chain consisting of N > > 1 mass particles which are assumed to be harmonically connected 
to their next neighbors. Any particle in the chain p (p = 0, .., N — 1) is assumed to have one degree of freedom 
characterized by the displacement u p (t), where t indicates the time coordinate. The Hamiltonian of this 
system can be written as 

N-l 
p=0 

The grading parameter £ characterize the exponential spatial dependence of the masses m p = mQ^ 2p and 
of the spring constants f p = m p uiQ. We assume that £ > 0, where £ = 1 represents the "trivial" case of a 
homogeneous linear chain. A schematic representation of the chain is drawn in Fig. 1. 

Our goal is to deduce the N x N frequency domain lattice Green's function containing the complete 
dynamic information such as the vibrational mode spectrum. The Hamiltonian equations of motion of the 
above system (1) are obtained from 

VHP u r = (2) 
where absence of external forces is assumed. This equation reads 

e p u P = -oji{e p {u P -u p+l )+e p ^\u p -u p ^)} , p =o,..,n-i (3) 

where strictly speaking this equation is not defined at the boundaries p = 0, N — 1 without imposing boundary 
conditions, i.e. assumptions about the fictive displacements u-±,un- In the case of a homogeneous chain 
(represented by £ = 1) it is convenient to impose periodic boundary conditions. Let us for a moment disregard 
boundary conditions to re-introduce them later in a quasi "natural" way. At present (3) can be rewritten 
more compactly in matrix form as 

ii = -A" 1 K u (4) 
Here we introduced the (dimensionless) N x N mass-matrix A 

A P9 = S pq f p (5) 
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where S pq denotes the Kronecker symbol (no sum over p is performed) and the N x N stiffness matrix K 
which is given by 

K Pq = e p "l (<W1 + ^) - Sp+i, q - (6) 

Whether or not K is positive definite or positive semi-definite, i.e. whether the eigenfrequency zero 
exists which refers to a uniform translation u p = at + b Vp of the chain, depends on whether or not such 
uniform motion is compatible with the boundary conditions. We have to note that the matrix A _1 K of (4) 
is non-symmetric, i.e. non Hermitean. Hence its eigenmodes are mutually non-orthogonal. 

To determine the eigenmodes and eigenvalues of the "dynamic matrix" A _1 K it is convenient to introduce 
an auxiliary vector y defined by 1 

y = A5 • u (7) 

having the Cartesian components y p = i p u p (no sum over p). The equation of motion for y assumes then the 
form 

y = -Ly, L = A^Ka4 (8) 

Matrix L has the components L pq = £~( p+q }K pq (no sum over p, q) with 

2 

L pq = ^ (Ml + £ 2 ) " aS P +i, q + Sp-i, q )) (9) 

We observe that L pq = L qp is a symmetric, Hermitean (self-adjoint) matrix which coincides in the de- 
generate case £ = 1 with K pq representing the homogeneous linear chain. Hence the Hamiltonian (1) can be 
rewritten as 

(N-l n-i \ 

E Vl + E LpiVpV, (10) 
p=0 p,q=0 J 

The rearrangement of terms which lead to expression (10) is justified when we assume periodic boundary 
conditions for y p , namely 

y P {t) = y P +N{t) (11) 

Let us further take into account that 

L = A^ (A^K) A-5 (12) 

From this relation follows that L has the same spectrum of eigenvalues as the dynamic matrix A _1 K which 
characterizes the motion of the "true" displacements u p = i~ p y p - Hence we can determine the frequency 
spectrum of A _1 K by just determining the spectrum of L. Due to the symmetry of L it has has mutually 
orthogonal eigenvectors (unlike A _1 K). 

Equation (9) suggests that the eigenvectors v of L have components of the form 

v p = Cfi p (13) 
with C being a constant and the admissible values of /x are yet to be determined. Then we have 

N-l 2 

E L pq v q = CV |f {(1 + £ 2 ) - £( M + /x- 1 )} = Xv p (14) 
From this relation we find that to these eigenvectors correspond obviously the eigenvalues 

a = f {a + a - eoi + p- 1 )} = f {« - aok - ^)} (15) 



1 By this substitution we transform the problem into a self-adjoint one with respect to the new variable y. 
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Taking into account the periodic boundary conditions (11) which must be reflected by the eigenvectors to 
fulfill these boundary condition, namely v p = v p+ n, leads to the condition for fi as 

» N = 1 (16) 
Hence the admissible \i are just the N complex unity roots given by 

fi m = e ikm , k m = ^^, m = 0,..,N-l (17) 

By imposing periodic boundary conditions (11) we define at the same time L in (9) at the boundaries for 
p = 0, N — 1 by its periodicity condition 

L pq = L p+Ns , p = 0, .., N - 1 (18) 

which holds for both indices p and q because of its symmetry L pq = L qp . With (11) both L and K become 
completely defined N x N matrices, i.e. for £ ^ 1 positive definite 2 and for £ = 1 positive semi-definite. The 
N normalized (Bloch) eigenvectors v*"^ (m = 0, ..,7V — 1) of L have the components 



(m) _ ikmP ciq\ 

v p - r^e w 



fulfilling the condition of ortho-normality 



N-l 



E = 5 ™n (20) 

p=0 



and the completeness condition 



7V-1 



E 4 m) v q m > = S pq (21) 

m=0 



With k m = jj-m the explicit N eigenvalues read 
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\ m = to 2 m = ^{(ti-cosk m ) 2 + sm 2 k m } , m = 0,..N-l (22) 

2 2 

with the dispersion relation uj m = uj{k m ). For £ / 1 it follows that < ^(1 - £) 2 < A m < ^(1 + £) 2 . The 
frequency ojq = occurs for &o = only in the degenerate case of a homogeneous chain when £ = 1. By 
introducing the periodic boundary condition (11) matrix L has for ( / no eigenvalue equal to zero. This 
is due to the fact that (11) does not admit uniform translations in all components u p = £,~ p y p which would 
correspond to an eigenvector equal to zero. The following further observations can be made: 

(I) Limiting case £ >> 1: This limiting case yields for all eigenfrequencies lim^oo uj„ l = loq where the 
band width tends to zero as 2^. By increasing the grading parameter one can "squeeze" the eigenfrequencies 

u) m all being localized within the band u>q(1 — |) < uj m < u>o(l + |). 

(II) Limiting case £ — > 0: lim^o^m = ^ — > oo where the frequencies are localized within the band 
u;o(| — 1) < u> m < o>o(| + 1) of constant width 2ooq. By decreasing the grading parameter we can "drive" the 
eigenfrequencies to arbitrarily high values where the band width remains constant 2ujq. 

The possibility to control the band of eigenfrequencies by the grading parameter may open interesting dynamic 
applications. From (12) it follows that A _1 K has the same eigenvalues (22) where its eigenvectors have the 



"Since the uniform translation u p (t) = at + b is not compatible with (11) 
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components £ v v p m ^ = ^f e * fcmP ( m = 0) ">-^ ~~ 1)- We observe that the p th components of each eigenvector 



v m are scaled with the same scaling factor £ p which is therefore true for any admissible displacement u p 
which thus takes the modal expansion 

p-v N ~ 1 

Up {t) = C p y P (t) = ^Y, ( Ane i(fcmP_Wmt) + S m e i ( fcmP+Wm *)) (23) 

Equation (23) represents the general solution of the homogeneous 3 system (4) with 2N arbitrary constant 
coefficients A m , B rn which can be uniquely fixed to satisfy any 2N initial data u p {t = 0),u p (t = 0). The motion 
associated with (23) is hence characterized by the Bloch waves y p (t) being spatially exponentially deformed 
by the factor £~ p . From (23) we observe the scaling relation u p +jv(i) = ^~ N u p (t). 

For the following deduction of the spectral oscillation mode density density it is convenient to employ the 
apparatus of dynamic Green's functions. Let us first define the time-domain Green's function Q(t) describing 
the motion of y(i). To this end we consider the inhomogeneous differential equation system 

{(| + 6) 2 1 + L}.y = b(t) (24) 

where e > denotes an infinitesimal positive damping constant, 1 = (8 pq ) the N x N unity matrix, and b(i) is 
an arbitrary source term which determines the response y(i) uniquely. Expressing the solution y(i) in terms 
of a convolution yields a representation of the form 

/oo 
G{t-t')-h{t')dt' (25) 
-oo 

where the kernel Q is the (time domain) Green's function of (24). By taking into account the trivial identity 

/oo 
5(t-t')l-b(t')dt' (26) 
-oo 

thus (24) requires that 

{l(! + e ) 2 + Ljc7(t) = 5(t)l (27) 
which defines the Green's function Q{t). By introducing the Fourier transforms 



2vr 

and 



Qit) = i- / g{u)e~ ibjt du (28) 

oo 



1 



'OO 



S(t) = — I e- iujt du (29) 
2ir J-oo 

we can transform (27) into the frequency domain to determine the dynamic Green's function Q{oj) in the 
frequency domain (which is nothing but the Fourier transform of the time domain Green's function Q{t)) to 
arrive at 

Q(uj) = j(L - l(u + ie) 2 | 1 (30) 

where Q pq = Q qp has the same symmetries as L and contains the complete dynamic information on the system, 
especially on the oscillation spectrum (density of normal oscillators) which we will deduce below. Note that 
we introduced an infinitesimal damping constant e > which guarantees that the Fourier representation (28) 
of the time-domain Green's function Q(u) remains well defined, i.e. L — 1 [u> + ie) 2 remains invertible along the 
entire real cj-axis, especially also at the eigenfrequencies co = uj m . In the following any expression of the form 



'Homogeneous" due to the absence of external forces. 
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f(u + ie) is conceived as its limiting case e — > 0+ where the zero is approached from the positive side e > 0. 
Inverting (27) requires such an infinitesimal damping constant which regularizes the problem and guarantees 
the time-domain Green's function to be causal [8] (i.e. Q(t) = @(t)Q(t) where Q(t) denotes the Heaviside 
unit step function). 

Before we evaluate (30) explicitly it is worthy to consider its relation to the (non-symmetric) Green's 
function of the "true" displacements u(t) which is defined by 



G(c 



A _1 K- l(u + ie) 



The transformation between G and Q is easily gained by rewriting (31) in the form 



G(w) = [a _ 3 KA-5 - i( w + ie) 2 ) A3 



or 



G(w)=A - 5 A" KA 2 - i( w + i e ) 2 



A2 



hence 



G(w) = A~5 5( W )A2 
G(cj) is a non-symmetric matrix having the components 



(31) 

(32) 

(33) 
(34) 

(35) 



where no sum over p, q is performed. The same relation connects also the Green's functions of the time 
domain. Let us evaluate now (30) in explicit form. To this end we represent this Green's function in its 
spectral representation 



N-l 



= E 



Am) 



m=0 W m 



(u> + ie) 



v 



(m)* 



(36) 



where * means complex conjugation and v^ m ) denote the Bloch-eigenvectors of L and (a<8>b) M = a p b q indicates 
dyadic multiplication. Hence 

N ~ 1 p ik m (p-q) 



Gpq(uj) 



1 

N 



E 



m=0 m 



uj1 1 -{u + ie) 2 



(37) 



Within the frequency band the Green's function (37) appears to be complex valued due to the presence of the 
imaginary infinitesimal quantity ie. This infinitesimal term becomes irrelevant outside the frequency band 
where the Green's function is purely real. In view of the symmetry of Q pq one can conclude that Q pq = Q\ p - q \- 
To evaluate (37) explicitly we use the relation 



i N ~ 1 i r 2n 

Y,f{e ik ™) = - f(e* k )dk 



(38) 



which holds asymptotically for N >> 1 where k m = ^p 1 G [0, 2ir] (m = 0, .., N — 1). Hence by using (38) we 
can write the Green's function (37) with the explicit dispersion relation uj m = co(k m ) of (22) as the integral 



tpq 



1 

27 



2tt 



Jk(p-q) 



uj 2 {k) - (uj + ie) 



-Ak 



To evaluate (39) explicitly by employing dispersion relation (22) we set 



u 2 {k) + ie) 2 = - e+){ji - e 



(39) 



(40) 



with n = e lk and 
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{n 2 D + n 2 -2u 2 ) 

a = cosh ^ = (n^-ng) (41) 

where we introduced the lowest eigenfrequency 4 $7o = — C\ an d the highest (the Debye-) eigenfrequency 
^(1 + £) = n D of the chain. 

The argument = x + iip i n this equation is generally a complex quantity. In the case that the zeros 
e ±<t> f (40) are located on the unit circle, we have to replace in (41) u by u + ie. This is the case when uj is 
within the frequency band. To further evaluate (39) we write it as a complex integral over the unit circle in 
the complex //-plane (where fj, = e lk , dk = (z/x) _1 d//) 

* 1 M -d/z (42) 



2™ 2 J H=1 (/i-e^- e -^ 

To evaluate this integral by using the theory of complex functions we have to determine the residuum of 
the integrand within the unit circle = 1. Without loss of generality there is always one argument <p = 
with x > such that the pole at is located within the unit circle since |e~^| = e~ x < 1. Hence evaluation 
of (42) yields Q pq = 2iriRes(..)\ lx=e -<t, to arrive at 



£ e -\p-q\<t> 

2uo 2 sinh 

To write the Green's function explicitly the following three cases have to be considered, namely 



5«=i-rr (43) 



Case (i): tt = f\l - £\ < u < f(l + Q = £l D 
Case (ii): < lo < $7o 
Case (iii): t-o > £Id 



Case (i): This case is relevant for the determination of the density of vibrational modes, i.e. when 
< w < is within the band of eigenfrequencies. In this case we observe that a of eq. (41) takes 
the values o,(Qd) = — 1 < a < a(Qo) = 1. Hence we can rewrite (41) as 



(n 2 +n 2 -2iv 2 , 

a= (n 2 D -n 2 ) = COS( ^ o<^<tt (44) 



since </?(f2o) = and (p(Qo) = 7T- In order to avoid the singularities of the integral (42) being located on the 
complex unit circle we have to replace in (44) co by uj + ie with an infinitesimally small positive e which is 
equivalent to replace (p by <p + ie in the limit e tending to zero. Integral (42) is then well defined, since the 
singularities at e ±i</3 are shifted infinitesimally onto e izi ^ (fi+ie \ By this regularization we have to put in (43) 
(j) = e — ip (\ = e ^ 0+) to arrive at 



% e i<p\p-i\ 
'"' ~2uM~sm(p 



or 



<*(0 - JL (£±iVTEZt![ (46) 
where we have used that sin ip = — a 2 > in the interval < 92 < 7r which assumes 



4 We emphasize that fi = 11 - £| > for 1. 
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sin ^ = (n^-ng) V ( ^ " w2)( ^ 2 " ng) (47) 



With (47) and (44) we can rewrite the Green's function (46) in the form 

= fc^ + ng - 2^ 2 + 2iy/(nl - ^){^ - n 

^-^)( w 2_ n 2) _ ng) ^ 



(48) 



where (48) is restricted to case (i), i.e. to frequencies within the frequency band Q,q < oj < Q,£>. In this case 
we also obtain the density of vibrational modes by using the relation [2, 4, 9] 

2uj 

p(u) = — Im TrG(uj + ie) (49) 

7T 

In (49) Tr(..) denotes the trace of a matrix. We obtain for (49) 

ttuJq sin <p ttuJq y/l - a 
to arrive finally at 

p(oj) = = , Slo < oj < Qd (51) 

which is defined and nonzero only within the frequency band and zero elsewhere since the Green's function 
outside the band is purely real (cases (ii), (iii) below). 

An important observation is that (51) possesses two singularities, one at the Debye frequency like 
in the case homogeneous linear chain, and a second one at the lowest frequency J)o which disappears in the 
degenerate case £ = 1 of a homogeneous chain. This second singularity at is due to the suppression of all 
modes with frequencies u> < Qo which are shifted to reappear at frequencies uj > CIq. 

The mode density (51) is normalized such that [2, 9] 

n D 

p{uj)duj = N (52) 

n 



To check this relation for (51) we have to confirm that 

<" n D 2Nlo 



f 



-.du = N (53) 



It is convenient to use relation (47) to evaluate 

/ ^)dw = 4 — ^T dw ( 54 ) 
Jn kuJq Jn Q smip{L0) 



To evaluate this integral we can directly make use of the substitution (47) together with (44) indicating 
that q z 4i ^qz duj = simpdip. Then we confirm the required relation (52), namely 



p(u)dio = -(ip(n D )-ip(n )) = N (55) 

where from (44) follows that (p(£lo) = arccos (1) = and <p(Qd) = arccos (— 1) = it. The mode density (51) 
degenerates in the case of a homogeneous linear chain represented by £ = 1 where Qo = and = ujd = 2cuo 
to 
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Pt=i(w) = — : , < uj < uj d (56) 

which is the well known classical expression for the oscillator mode density of the homogeneous linear chain 
which can be found in any textbook of lattice dynamics (e.g. [2]). Expression (56) can be also reobtained 
directly from the dispersion relation of the linear chain using p^ = ±(uj) = 2|^|^r;| and (22) 5 for £ = 1. 



Case (ii): 

In the range of frequencies < uj < Qq, i.e. outside the band of eigenfrequencies the coefficient (41) takes 
values a > 1 thus we can put 

(n 2 n + - 2uj 2 ) 

a= °4-"S) coshx ' x>0 (57> 

where <p = \ > is purely real. Hence the Green's function (43) assumes in this case 

£ e -\p-q\x 

Qf) = % c . (58) 

2uj\ sinhx 



which is purely real valued. In this expression sinhx = Va 2 — 1 is the positive root because of x > 0. We 
thus have 



sinhx = (^-n§) V ( ^ " w2)( ^ " w2) (59) 

Expression (58) hence has to be rewritten as 

= _ L= ^ ^ + ^-^ :! y^- r ^-^ b - g| (6o) 



(n| + n 2 - 2u? - 2y/(a* D - gx^j - ^ 2 )) 



where < uj < We note that (60) contains the static limit Qpq (uj = 0) = (L 1 ) pq which is obtained by 
putting uj = 0. 

Case (iii): It remains to discuss the case when the frequency exceeds the Debye frequency uj > fljy. In 
this case we have 

(Q 2 n + nZ - 2uj 2 ) , 

° = (n^-ng) - = cosh(x + t7r)<-i, x>o (6i) 

that is a = — coshx- Hence we have to put <p = x + in in (43) to obtain the Green's function as 

...... (_l)p-q+l£ e -\p~q\x 

y ™ ~ 2uj 2 sinh X { ' 



Since x > the sinhx has to be as in case (ii) the positive root, namely sinhx = \/a 2 — \. Thus we have as 
in case (ii) 



smhx = {n 2 D _ Q 2 ) y/^ i D-^W-^) 



(63) 



5 where the prefactor "2" accounts for the number of branches of the dispersion relation ui(k) in the fe-space. 
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To rewrite (62) we have to note that coshx = —a where a is defined by (61). Taking this into account 
the Green's function has to be rewritten as 



= (-D (ng, + n§- ^ + 2^ -^) ( n§-^)) |p ~ g| 

where u > 



3 The Continuum Limit 

We consider in this section the transition of system (1) to a continuum. To this end we assume all particles 
having the same distance h to each other with h — ► . We assume the chain has length L = Nh where L 
is kept finite and constant in this limiting process. As a consequence the total number of particles becomes 
infinite N = ^ — > oo. Particle p has then spatial coordinate x = ph where p = 0, ..N — 1 varies x quasi- 
continuous. The displacement u p (i) of particle p has then to be replaced by the displacement field u(x,t) 
where the the displacement of the "neighbor" particle p + 1 is given by u(x + h,t). In the continuum limit 
any sum over p can be replaced by an integral, namely 

-'■ Ax 



EF,~l (65) 



p=0 

where L = Nh and dx « h. This relation becomes asymptotically exact for a small h — > 0. We further 
introduce [11] 

^p+i(i) — w p (t) = u(x + /i) — u(x) ~ h—(x, t) (66) 

where we take into account only the linear order with respect to h. The Hamiltonian (1) can then be rewritten 
in this limiting case as the integral 

*-3jG l *i(£) , +*«(5)> 

where we have set £ 2p = e 2/3x with p = | and /3 = lim^^o ^jf- This is justified because we have to demand 
continuity of the function £ p ~ £ p+1 , i.e. ^ ~ 1 when h — ► 0. Hence it is justified to set in this limiting case 
£ = e' 3 ' 1 1 + (3h where j3 is constant. In this way we can obtain the limit of £ p = with p = | for a given x 
indeed in the form 

x 

i p = \mx{l+f3h)h =e Px (68) 



Since the integrand of (67) represents a finite energy density we have to demand in the limit h — > the 
relations 

mo(/i) = Poh , wo W = ^2 ' P = ~j~ ( 69 ) 

where we introduced the constants Q, po and (3. Then the energy density (energy per unit length) in (67) can 
be written in the form 

n = P M^(duy + ^(duy (?o) 



dtj 2 \dx) 

The integrand represents the Hamiltonian density where the mass density pm{x) and elastic modulus 
p(x) = Q 2 pm(x) are exponentially graded 

Pm(x) = po e 2/3x , p{x) = n 2 Po e 2l3x (71) 



10 



By conceiving integral (67) as the energy functional, we obtain the equation of motion for the unknown 
field u(x, t) by employing the Hamiltonian principle [10] for fields, i.e. by putting the functional derivative 
^ = to arrive at the partial differential equation 

^S=^(>^ (72) 



dt 2 dx V dx J 

which corresponds to the case of absence of external forces. This equation can be also obtained directly from 
(3) by performing the transition by using (66). Hence (72) is the continuous counterpart of (3). In order to 
resolve (72) we make the transformation 

u(x,t) = y{x,t)e- px (73) 
which is the continuous version of (7). By this substitution (72) takes the form 

d 2 y 

^ = -Cy(x,t) (74) 



where we introduced the differential operator 



Equations (74) together with (75) can be written in the form 

y + P 2 y = o (76) 



Q 2 dt 2 dx 2 



which has the form of a Klein- Gordon equation of the scalar field of a free moving (spin-0) particle [12]. 
Further we introduce the periodic boundary conditions 

y(x,t)=y(x + L,t) (77) 

which are the continuous counterpart to (11). 

We show now that (74) is the continuous counterpart of (8) where the differential operator 
C is nothing but the continuous counterpart of the matrix L of (9). To this end we consider 



E^" 1 L pq y q (t) = lo 2 {(1 + e-^)y p (t) - e^ h (y p ^(t) + y p+1 (t))} 

= u 2 {(1 - e- h ?) 2 y p {t) - e-? h (yp-i(t) + %>+i(i) - 2y p )} 



(78) 



Taking into account only terms up to order h the first term in the last equation contributes with P h uj^y(x, t) = 
£l 2 (3 2 y(x,t) and the second term with —h 2 L0 2 -^y(x,t) = —Q 2 -^ 7 y(x,t) to arrive finally at 

N-l o 

E l pivS) ~ ~ n2 (g^ - ^M*- *) = £ y(*< *) ( 79 ) 

The normalized eigenfunctions of operator C can be written as the continuous version of the Bloch eigen- 
vectors (19) 

W M( S ) = _L e ^-* j K m = (80) 

with m = 0, ±1, ±2, .., ±oo. The eigenfunctions v^ m \x) fulfil the periodic boundary conditions (77). The 
dispersion relation is then obtained by the eigenvalues of C, namely 

Cv {m \x) = Lo 2 {K m )v^ rn \x) (81) 
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and yields 



u; 2 (K m ) = n 2 (K 2 m + (3 2 ) 



(82) 



This dispersion relation has to coincide with the dispersion relation (22) in the continuum limit. To this 
end we rewrite (22) where we note the relation between the wave numbers k m = hK m and obtain 

u4 = u l(e iK ™ h -P h - i)( e -i K mh-ph _ ^ w u %h 2 {Kl + (5 2 ) = n 2 (Kl + (5 2 ) (83) 

which recovers (82). We note that the dispersion relation (83) represents a discrete set of eigenfrequencies 
due to the discrete values of K m . This is a consequence of the periodic boundary conditions (77) with a finite 
length L of the chain. In the limiting case L — > oo the spectrum of eigenfrequencies becomes quasi-continuous. 



Continuum limit of the Green's function. In the following we deduce the continuous counterpart of 

Green's function (43) and above cases (i)-(iii). Then we have = ^r^o = ^o(e l3h — l)e _/3/l ps [3ujq}i = (3Q 

and in the same way £Id ~ 2^ + [50. — > oo. This observation, that the lowest frequency is at £Iq = with 
no finite Debye frequency is also reflected by the dispersion relation (83). We have hence in the continuum 
limit only the two cases: 



Case (i): f\) = (3£l < uj < oo 
Case (ii): < u < Q = 0£l 

To consider these two cases it is useful to write coefficient (41) in this approximation 

« = l + y(/3 2 -^) (84) 

where the lowest non- vanishing order in h is the second order term which has been taken into account. In 
view of (84) it follows that in case (i) a < 1 and in case (ii) a > 1. 

Case (ii): Here we can put 

a = CO sh X = l + y(/? 2 -^) > 1 X >0 (85) 

since P 2 — > and by taking into account only the dominating order in h we have 



s inh X = /*yV-^J (86) 

where the positive root has to be taken since x > 0- Let us consider the expression (58) for Qpq in the 
continuous approximation. Here we have to put p — q = xx ~ h £2 = f and the fact that a matrix multiplication 
with the Green's function becomes a convolution integral where the matrix of the Green's function becomes 
a convolution kernel. Since the summation of a matrix multiplication has to be transfered into an integral 
by taking into account (65), we have to multiply the discrete Green's function with a factor j- to arrive at 

gd^l, lo) = lim^o {j^Qpq^j ■ Then we can generate the continuous limit of the Green's function by constructing 
e~ x from (85) and (86) and by using the Green's function matrix of (58) in the form 

\x\ 
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which simplifies when we take into account the relation 



lim (l + Cl h + c 2 /i 2 ) h = e ClX (88) 

where the quadratic term in C2/1 2 and higher powers in h are irrelevant. Utilizing this relation yields the 
Green's function in the form 



1 e- W V^ 

g (M) (N^) = ^o2 / (89) 

where < uj < (50. It is important to note that because of the periodic boundary conditions (77) which are 
to be fulfilled, the Green's function has the form (89) in the interval < \x\ < L and has to be L-periodically 
continued. In the limiting case L — > 00 (89) is defined for all < \x\ < 00. 

Case (i): In a similar way we generate the continuous limit in the band of eigenfrequencies, i.e. when 
(50 < uj < co. In this case we have a < 1 and put 



h 2 



UJ 



2 



a = cos^ = l-'—q--^ 2 ) < 1, ip>0 (90) 



and 



sin^ = W— -(5 2 (91) 



where the positive root has to be taken since tp > 0. In a complete analogous way as above in (87) we consider 
the limit h — > of (45) by putting g = \Qpq and obtain 



1 e 



,2 



9 {i) (\A^) = —t ~ a (92) 

where co > (30.. By utilizing (92) we can also generate the continuous limit of the mode density which is non- 
zero within the band of eigenfrequencies u > (50. To this end we have to employ the continuous counterpart 
of (49) to arrive at 



which yields 



p{u) = —Im g^(\x\ = 0, to) (93) 

p(u) = 1 ^== (94) 

As in the discrete case this mode density is only defined within the band u> > (50 and zero outside for 
< uj < (50. Expression (94) is also in accordance with the continuous limit of (51), by using the asymptotic 
relations f^o = (50 and Op = 

The deduced Green's function g(|a;|,u;) of equations (89) and (92) is the fundamental solution of the 
frequency domain representation of (74) and fulfills the equation 

(C-uj 2 )g(\x\,uj) = 6{x) (95) 

where 5(x) denotes the Dirac (5-function. This equation is noting but the continuous analogue of the matrix 
equation leading to (30) in the discrete case. We note that when we impose periodic boundary conditions the 
Green's function (89) and (92) is defined at < \x\ < L and has to be //-periodically continued. In the case 
L — ► 00 of an infinitely long chain these expressions hold for the entire range < \x\ < 00 and represent the 
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Green's function of the ID infinite space. Equation (95) can be rewritten in the form of an inhomogeneous 
Helmholtz equation, namely 

-tf(Jj + £-^9(\*M=SW (96) 

Obviously the character of solution of this equation depends on the sign of ^ — (3 2 reflecting the two cases 
(i) and (ii). 

It remains to verify that expressions (89) and (92) indeed solve (96). To this end we consider 

rl 2 

-f- ce A N = A 2 c e A|a;| + 2c A Six) (97) 
dx z 

where 0(x) denotes the Heaviside unit step function and c denotes a constant. We further have used 

^ = 1 {@{X) ~ @{ - x)) = 26{X) (98) 
and 5(x)e x \ x \ = 5(x). We hence can rewrite (97) in the form 

(^ - A 2 )ce AN = 2cX5(x) (99) 



32 w 

TP 



By putting c = jxn? we recover with A = i\j ^ — ~jp expression (92) of case (i) and with A = —\ff^ 2 
expression (89) of case (ii). We note that solutions with A's of opposite signs also solve (96), however they 
correspond to physically inadmissible 6 solutions of (96). The continuous counterpart to (35), i.e. the Green's 
function of the "true" displacement field u(x,t) is given by 

g{\x\,u)=e-f >x g(\x\,«;) (100) 

where g(|x|,u;) is the Green's function given by (89) and (92), respectively. We further observe that (100) 
fulfills 



- n ^r 't) + ^ s } =s ^ (ioi> 

which is the continuous counterpart of the matrix equation defining the discrete Green's function (31). 



4 Conclusions 

Expressions (48), (60), (64), respectively, provide the explicit representations for the dynamic N x N Green's 
function defined by (30) for the entire range of frequencies. This Green's function is linked with the Green's 
function of the "true" displacements simply by (35). The degenerate case £ = 1 recovers the well known 
expressions for the Green's function of the homogeneous linear chain which have been known for a long time 
(see e.g. [2, 9]). 

We also deduced the continuous limit of this system which corresponds to a ID graded linear elastic 
material. The transformed equations of motion take in the continuous limit the form of a Klein-Gordon 
equation (76) where the limiting transformation of the Green's function yields the expressions (89) and (92), 
respectively, which solve (96). 

In spite of the simplicity of the model system, the exponentially graded linear chain and its continuous 
counterpart seem to be one of the simplest systems representing one-dimensionally graded materials which 
are accessible to a concise analysis. 

6 A = —i\J^z — /3 2 leads to a non-causal Green's function in the time domain in case (i) and A = \J P 2 — to a divergent 
solution for L — > oo at x — > oo in case (ii). 
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The band of natural frequencies, i.e. its width as well as the highest and lowest frequencies can be 
controlled simply by the grading parameter £. This suggests the possibility to design exponentially graded 
FGM with frequency bands according to desired dynamic properties. 

Possible further extensions of the present analysis could be on linear chains with self-similar (fractal) 
properties [13]. We hope that the present analysis inspires also work in such directions. 
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